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XIII-3. 


Dopo lo studio degli operatori Fuchsiani iperbolici fatto da Hide- 
toshi Tahara [T], Bove-Lewis-Parenti [BLP1-2], e altri, vorremmo considerare 
la soluzione numerica mediante elementi finiti di un sistema iperbolico Fuch- 


siano: 
(+) Lu(t,x) = ta.u + tA(t,x,3, )u+B(t,x)u = f(t.,x) 
ove supponiamo che 


(i) L'operatore aptA(t,x,9,) sia un sistema NxN di operatori dif- 


ferenziali, simmetrico e strettamente iperbolico; valga la disuguaglianza 


|Re(Au,u)| do luj?, 


N 
UECS(RI). [(u,v)= {Patata = (u,u)]. 
R" 1 
x 


(ii) B(t,x) sia una matrice di funzioni che soddisfi 
Re(Bu,u) ® b pu]? b >O. 
? b) ? o 


Sia feco (RARE), ci sono soluzioni singolari di (+), ma c'è una 
soluzione sola che sta in C°(10,T]; 2'(R7)); inoltre il dominio di dipendenza 
della soluzione liscia di (+) è il solito "cono di luce". Se f ha una singola- 
rità C°, l'insieme WF(u) è contenuto (per t=0) nelle strisce bicaratteristi- 
che. Per dettagli, facciamo riferimento a [BLP1] e [T]. 

Presento, seguendo Tahara, qualche stima dell'energia per (+), e 
poi presento due metodi di approssimazioni numeriche delle soluzioni di (+). 


La strada da seguire è la seguente: 
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(1) Porre (+) in una forma debole e fare stime di energia; 
(2) Introdurre la forma debole su uno spazio di elementi finiti nello spazio 
e cercare di fare il confronto dell'errore con l'errore di una certa proie 
zione della soluzione; 
(3) Fare una discretizzazione anche nel tempo e dare una approssimazione "pie 
namente discreta" della soluzione. 
In un recente lavoro, Anatoly Genis [G] ha portato avanti questo 


programma per l'operatore di tipo Euler-Poisson-Darboux: 


Lt) = talu + (2p+1)au 21008, (8,008, 1)-c lau) = 0 


u(0,x) = g(x). 


Genis ha ottenuto risultati anche di stime in L°, superconvergenza, 
ecc. Il nostro studio è ispirato a Genis. I metodi adottati qui sono modifica- 
zioni dei metodi di Todd Du Pont [Du2] e Lars Wahlbin [W]. Essi (nel caso non 


Fuchsiano) hanno risultati anche per equazioni quasilineari. 
STIME DI ENERGIA PER IL PROBLEMA CONTINUO 
Ci interessano risultati locali nello spazio; per semplicità la- 


. o psé n _ 
voriamo nel caso periodico su Tr Per u = (Use aUn)» v= (AETEEIAZO, deno- 


tiamo con (u,v) il prodotto scalare in 


2 N ù 
L (T): (u,v) =)_ J uv da. 


1 7 
x 
2 
Per r=0,1,2,..., |u|.= È Jada dx; 
la|=r 


La 
= » Lx 2. = 
Du} = Duta 5 dalf =D Tal 3 Jul, = sup lu0o1. 


Su G consideriamo il problema 


(+) Lu = ta u + tA(t,x,3 Ju + Bu = f(t,x) 


con L che soddisfa (i) e (ii). 
Poniamo (+) nella forma debole: 


(+), (taju + tAu + Bu, v) = (f,v), 


Vv, per esempio, in cleto,ti;#t (7). 


Con v = u(t): 


1 
i 


34 2 u(t)1? - a.tlu(t)1 + b_Ju(t)1? s 1F(t)] lu(t)]. 


Dalla disuguaglianza ab s & i? *% È bi, risulta 


48 


2 


tdi lu(e) 1° + (6 2a t-0)lu(t) 1° s T 1e(0)1?. 


le 


Fissiamo T piccolo tale che = bd = bo - a 7990. 


mie 


Per ogni t < T, 


- - 1 
-b, d ,,b 2. 2 
total) sf 


Segue 
t=t 


= t - 
t°lu(t) |? | scf sP"1/#(5) |Pas. 


t=e € 
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Sotto l'ipotesi lim tP|u(t)|f=0 (che taglia fuori la "soluzione singolare"), 


t+0t 


XIII-6. 
t Di 
Ju(t)1° = cf di UO] 
(0) 


Questa maggiorazione, sfruttando la disuguaglianza di Hardy, dà varie stime con 


peso # in LP(dt). In particolare sfruttiamo solo l'ipotesi B>0 per ottenere 


2 2 
(*) sup Jult)| sC sup |f(t)[C. 
OstsT OstsT 
La dimostrazione esposta dà (*) per T piccolo (sT,). ma per 
TsTs i » tale stima risulta dalle solite sitme per sistemi iperbolici. 


Passiamo alle stime per lu(t)], Prendiamo 8, dell'equazione (+). 
Risulta ì 


ta Ue. + pitti + ui derit + Bu, + ea = fe 


Sia V= (vu) = (au, sets, u), V soddisfa il "gran sistema": 
i n 


ta. V + tA(t,x,9,)V + Bv + tlv = ve + Du, 


t 
Ai 0 B! 01 
24-\ O, (-- i n 
dove È = QI AL-___k > 8= f 0! 8’ 
N ed eci Sl silla mcr 
LILLA. Me pesi di 


e È e D sono matrici (n° x n°) di funzioni. 


Anche se ovviamente DR non è strettamente iperbolico, si può fa- 


re la stima dell'energia di cui sopra; Du è già controllato e 
2 
Re((B + tC)v,v) = (bc) IVI". 


Otteniamo 


(*) 


2 2 
L4 sup|u(t)| ps € sup If(t)l, 


OstsT 
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La stima Pi serve quando facciamo un'approssimazione della solu- 
zione nello spazio. 
Per controllare apu(t), facciamo una derivata rispetto t dell'equa 


zione (+). Risulta 


(ta, + tA)u, + (B+1)u, + {tALu+AU+B,u} = fi 


Cioè: Sia w = du 


= - + a 
[ta, + tA + (B+I)]w fi {tAu + Au BU} 
Il ragionamento già fatto dà la stima 


sup |a,u| sc supt|f.J + Ju ,,} s C suptlf.]. + |fl.,,3. 
O<t<T ES Lai rtl tr rtl 


Riassumiamo le stime dell'energia ottenuta nel seguente 


Teorema 1. supla"ul s C supi] a6rl + nibti ..+[fl 
‘a CP EP t 


lt = rekl- 


SPAZI DI ELEMENTI FINITI 


Per una descrizione degli spazi di elementi finiti si veda Fair- 
weather [F] e Strang e Fix [SF]. 

Sia 2 = Lo (o più generalmente un aperto di R"); si considera un 
sottospazio gl di LE (0) di dimensione finita. Diciamo che sl è di classe 
MI, Osksr se e solo se 


I. she 4410) e 


II. Per 2sssr, per ogni ue H" (a), 
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inf (Ju-x|+h|u-x|,) s Ch°|ul_. 
e$ S Ss 
cab 
Di solito tale spazio % nasce da una sottodivisione di Qin trian 


goli o quadrilateri Te j=1,...,N,. Si suppone che ogni "elemento" ti ha dia- 


h° 
metro x h ed è regolare (cioè non troppo stretto, ecc.). Ogni vet ristretto 
ad un elemento Hi è un polinomio di grado sr-1 (o almeno un prodotto tensoria 
le di tali polinomi). 


La condizione sc peri 


(9) impone condizione di compatibilità sulle 
derivate di ordine sk lungo i bordi di due elementi con una faccia a vertice 
in comune. Nel caso n=1, la sottodivisione è semplicemente una partizione del- 
l'intervallo e x€ pet se e solo se xe ct. 

Per dettagli e esempi si veda la bibliografia. Citiamo esempi im- 


portanti: 


n=1: i polinomi cubici di Hermite (r=4, k=1); 


n=2: Su un triangolo, l'elemento cubico di Zlamal. 


Ad ogni vertice si dà XoXyXy? per trovare una base locale di po- 
linomi cubici (dim = 10),ci rimane un "grado di libertà", di solito si dà x 
al baricentro del triangolo. L'elemento di Zlamal è di classe MO: 
n=2: Su un rettangolo con lati paralleli agli assi si considera un prodotto 


tensoriale di polinomi cubici di Hermite. Per base locale (di dim = 16) su 


ogni elemento si dà LO MONO ad ogni modo per ottenere un sottospazio di 
classe Mi Il controllo di Xxy ad ogni modo assicura che le prime derivate 


normali sono continue attraverso un bordo comune di due elementi: cioè 


XE He (9). 
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IL METODO DI DU PONT NEL CASO SEMIDISCRETO 


Definiamo l'approssimazione semidiscreta di (+) la funzione 


U(t)e c°([0,T1;") che soddisfa: 


Per ogni VE si, 


(#), 
(ta U+tAU + BU,V) = (f,V). 


Attualmente la soluzione di (+), è la soluzione di un gran sistema di equa- 


zioni differenziali ordinarie. 


Notiamo che U(0) è univocamente determinato dalla relazione 
(B(0,x)U(0),V) = (f,v), ves. 

Inoltre 
b_|u(0)-U(0) |“sRe(8(u(0)-U(0)), u(0)-U(0)) 


= Re(B(u(0)-U(0)), u(0)-x) per ogni yes", 


dato che (B(u(0)-U(0)),V) = (f(0),vV)-(F(0),v)=0. 


Risulta che 
lu(0)-u(0)|sch"|u(0)],.- 
La stima per sup|u(t)-U(t)| è data dal teorema seguente: 
t 


Teorema 2. Supponiamo che s è di classe Lo e 


A = se lu(t)], + Itu,l, < 0, 


XIII-10. 
con u soluzione di (+) e U(t)e sh soluzione di (+), Allora 
Ju(t) - u(t)] < caten)n!"?, 
Prova. Sia W(t):[0,T] + sh un'approssimazione di u(t) che soddisfa 
lu(t)-W(t)] + hlu(t)-W(t)]sch"lu(t)], 
Jup(t)-W,(t)1 s ch']u,(t)],- 


: a NEEP sladigg i ; 
Per esempio, sia W la proiezione in H° di u: W è l'unica funzione 


in sh tale che 


g(U-H,V) = (u-W,V) + (v(u-W),wv) = 0, vesì. 


(Formalmente, W risolve approssimativamente il problema 
AW - W® Au-u. 


Tale problema è l'esempio fondamentale degli elementi finiti. Dato 8 non dipen- 
dente da t risulta che LA è la proiezione di u)- 
Scriviamo u(t)-U(t) = (u(t)-W(t)) + (W(f)-U(t)) = n(t) + c(t). Per 


ogni VE . 
(te, ttAt#Bn,V) = -(tn,+ tAn + Bn,V) s ca(n'+th""3)|v|. 


Sia V= c(t)e sh. 


Li Linteit” a (b_-a t-8)lc(t)1% sE a°cni” + 


C Ent” + pr? 
2° di 6 


h }. 


Segue, come nel caso continuo 


XIII-11. 


sup |e(t)|Zs cafen?"7(n%47%)). 


OstsT 


Finalmente |u(t) -U(t)|s|n(t)|+|c(t)]. 


Si chiede perché abbiamo preso una potenza di h. DuPont [Dul] ha 


dimostrato che per l'equazione 


u(0) = u(*) 
con i polinomi cubici di Hermite (qui n=1 e Ca di classe Nm. in generale non è 
possibile aspettarsi convergenza migliore di o(h°). 
DISCRETIZZAZIONE NEL TEMPO 
La relazione 
(4), (ta, U + tAU + BU,V) = (fV), ves! 
conduce ad un sistema (singolare!) di grande dimensione, il quale dev'essere 


risolto mediante un metodo di differenze finite in t. 


Sia F(t) una funzione (anche di x), definiamo, con At fissato: 


A = F(mat) ; Fn+1/2 = F((m+1/2)At); 
preze imm. se = mim 
2 ? tm At 


Usando la formula di Taylor con resto, 


2 
__m+1/2 (At) 1 
Un+1/2 UT 214 f (1- le})u,, {mt Atto > ? t) de; 
dat)? 1 st 
9 1Ume1/2 7 Sim © 314 Us |e1? + 3)st +0 3 )de; 


Per brevità scriviamo queste formule come 


m+1/2 _ (1)? 


Um+1/2 7 KUtt,m+1/2 > 


2 


su - (At) KaUttt,m+1/2" 


9tUme1/2 7 Sem 


Notiamo che 


si :6 sup |tv| , i=1,2. 


litoerl 


ltm+1/2 Ki Vm+1/2! 


Quest'ultima osservazione è di Genis. 
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Definiamo la soluzione approssimante "pienamente discreta" come 


una successione {U pP_ Cc sh che soddisfa 
m m=0 


pu(o) - u(o)] < ca n°! 
(4) nat Per m=0,...,M-1 


fee uni ve, 


m+1/2 m+1/2 h 
Uta 1/ 204 tte 28 me172) VO Bmg 3 Fnggza!) VES 
Sia 
3 m+1/2 m+1/2 
Lim = tre1/284m + tinta lm, Bmg» 


risulta che 


_ 2 
(LU) 41/2 7 Patt m (08) 9mt1/2? 


Im+1/23 Umeazztalttt,me1/2Ì tme1/20me1/2%10%t,m+1/2 


+ Baris2 bi'tt,mi/2” 


Possiamo portare gli operatori A sotto gli integrali (in t) che 


m+1/2° Èm+1/2 


definiscono Ki Risulta che 


lOme1/a! 5 SUP Mosa; + Itu, ly + lux 
m° mt1 


Con questi preliminari siamo in grado di enunciare una stima per 


Teorema 3.1. Supponiamo che Une sh sia già determinata e At,At/h sia 
no abbastanza piccoli. Allora la relazione 


_ h 
(Len?!) a” (fn+1,2?") 9, VEeEs È) 


determina U univocamente. 
mtl 


2) Supponiamo che 
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A = sup{|u], + |u,|, + |tu,,,|+|tu 
t 


ei eteltitu.zt +10, 1? 


Allora 
ta r-1 2 
un! s CA(1+10g(7)){(h+t)h + (At) "} 


Prima di dare la dimostrazione, notiamo che la costanie A può essere maggiora- 
ta del dato f. La discretizzazione nel tempo (almeno con questa dimostrazione) 
ci fa perdere un log(i/at). 


Prova. La determinazione di Und consiste nella risoluzione del 


1 
problema 


Per ogni ves” 


_ bt 
Brr+1/2Um+1/2?)=(tre1/2Umea*tme1/2 2 Amt1/2Ume1 


st = 
u 2 Br+1/2Un+1?!) ComoV) 


con p_ dato. 
m 


Per t+1/2 abbastanza piccolo 

Re B,101/2(VM)=ltna1/201- Sad, 2 vj? 
Se non c'è a priori controîlo su bo per tm+1/2 grande, sfruttiamo il fatto che 
Re B,43j2(V-V) > ile Sp +b, DI ivif. 


L'esistenza e unicità di Uni seguono dal Lemma di Lax-Milgram. 


1 
Per farci arrivare ad una stima per Una? prendiamo di nuovo 


sui as & ili ; 
W_ € si la proiezione in H di u . &W è la proiezione di Su . 
m m tum im 
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Sia Le a 


i r r 
ingl + ting # clunl * < CAh 


r r 
162%! S cis ih sCAN. 
Dalle relazioni 


2 
(Lem?!) = (f -(At) g 


m+1/2 m+1,2°") 


(Lit Un!) I EYZIAZE 


sia Und = (u-W)+ (MU) = Nn + En 


_ 2 
(L,#5ne") ni (tm + (at) Im+1/22")- 


Si nota che 


r r- 


1 Pi 
|L [= C*% + CAt_41/2"- + Ch A; l9n+1/2! s CA. 


atm m+1/28" 


Mettiamo V = mti/2e sì, 


t 


Ri m+1/2 
mt1/2 2 


2 t 
Stlemi + mensa lAmen/a8 , gMti/2, 


m+1/2 m+1/2 


+ (Bn+1/25 ,% ) 


s CAT. h° lan” 3 


m+1/2 
m+1/2 I 


2, _mt1/2 
+ ca(at)?|x9t1/2| 


miei 
[i 


Con3 = b - a t -6>0, risulta 
D) (ORaTI 
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:2 rr |_mt1/2,2 VAT 2,,2r-2 4 
tn+1/28+1%m! + bc | s CA {(tr+1/2)! )h +(At)} 
M-1 
Dividiamo per tn+1/2 = (m + 1/2)At e facciamo par At. Risulta 
M-1 M-1°, mt1/2,2 
pa s.1e,1%4t + b 3 le © [at 
m=0 mzo  (m+1/2)4t 
2,,2r-2 3° im hÉ ve A 
<s CA {h { t bt + } + (At) } 
+ 
7 mI/2 Tm) E (m+1/2) 
s cA°CtI h°"72 + (14109 M(h2" + (at)*)} 
bs 
La dimostrazione è fatta sotto l'ipotesi aotm < b78» per ty a , Îîl solito 
(o) 


Lemma di Gronwal] è applicabile. 


UN METODO DI WAHLBIN 


Per migliorare la convergenza in h, Wah]bin [W] ha proposto di in 
trodurre una "viscosità artificiale". Da ora in poî supponiamo n=1 e prendia- 
mo o lo spazio di polinomi cubici di Hermite su una partizione quasiuniforme 
di (0,1) (oT1). Ss" 


versa", cioè per x ES , 


è di classe Mi; in più soddisfa una cosidetta "ipotesi in 


hlxl,43 5 Clxlg» 4 3 051,23: 


g+l 


(Segue dall'omogeneità). 


A ? sa sue da 
Una conseguenza è una buona stima della proiezione in L di una 


. MET. 
funzione in H:: 
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Sia WES" tale che (u,V) = (W.V), ve 9, allora fu-w],sch*5]u],, 


s =-1;0,1,2. 


(Jul_; = SUP) a 


sel lol; 


Si veda Douglas-DuPont-Wahbin [DDW]. 
Per b fra 3 e 4 da scegliersi in un modo opportuno, sia u soddi- 


sfacente Lu = f, allora 


b b 
(tu, + tAu + Bu,V) # th fu, Va) (f.V) + th (u,,oV a 


xXx 
Definiamo U(t) € c°(10,T1;s") la soluzione di 


b 


(+) (tu, + tAU + BU,V) + th 


h,b 

Formalmente stiamo approssimando la soluzione di 
b 

tu, + tAu + Bu + th u = f; 


t KXXX 


un problema singolare di tipo parabolico. 
Teorema 4. Supponiamo che bo > 0 eb = 10/3. Mettiamo 
A = sup(|ul, + {tu,|,) 
t 
Allora 


sup |u(t)-U(t)]| < ca(n®+V/kp) 


ostsT 


2 


Prova. Prendiamo W(t) la proiezione di u(t) su sh in L. 


Scriviamo 
u-U = (u-W) + (W-U) = n+z 
Risulta (con V= x) 


) 


b 
(te, + tAx + By,c) + th (eno txx 


G b b 
= (tn, + tAn + Bn,c) - th (my Sax) + th (ui) * 


La parte sinistra è (come al solito) = 
ta 2 IRE 2 
2 dt SC) + (boat) ele) +th Je, 


Controlliamo la parte destra: 


Itnyo)] s s|z|? + ca°h; 


[OMolE s|z1Z + ca?h; 


b b 2 b,4,2, 
| th (mjet) Ss & th LA + Cth h'A; 


2 2b,2 
LthPfu, at) = IthP(u, ae) = stili? + ct nda, 


XX XX XXKX 


Il bello è il controllo di (tAn,z). 


I(tAnsc)]s Ctinl_,lcl, s ct Inl_,(l51, + ce DD) 


ctndacn"P/2pb/2 


A 


la, + 30) 


10-b 2.102 


n° 46 le + CENA 


US 


sth|z15 + Cth 
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Risulta 


t di i_2, di 4,2 b 2 
paglel +75 + (1-36)th al 


< calenBetn**Pientdatn10-0, 4219, 


Se b=10/3 , 10-b= 2b e 
suple(t)| s calen8 + en). 


Nota. Sfruttando una proiezione diversa sì può migliorare il ri- 
sultato fino a b = 7/2. Si veda Wahlbin [W]. 

Adesso seguiamo il metodo di viscosità artificiale di Wahlbin nel 
caso "pienamente discreto". 


Si considera una successione Wie es', Un? u(mAt), che soddisfa 
b 
lu(0) - ul s CA h 


Per ogni VE s 


(+) 
h,b,At 

dii b, m+1/2 
CL, Un?!) A ta1s2? ( x Kage 


VALE 


Teorema 5. Per At, At/h abbastanza piccoli la successione {Un} 


che soddisfa (+) è univocamente determinata da UL: Sia 


h,b,At 


= + 
A suprlul, Itu,ly + Itu,.l; + Itu,z[+lu}1}» 


allora 
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S LA 
* RO tm 178,4 2d 
lu Un! s CA {[1+log( a )1(h+(At)")+h} 


Prova. L'esistenza di U, segue come di sopra dal Lemma di Lax-Mil 
gram. 


. I ini aste ; 
Sia LIA la proiezione in L' di Un» Spezziamo 


usfu -W)+(W-u)3=n +%. ConL 


PA Le la n At * In+1/2 definiti come nella dimostra- 


zione di Teorema 3, 


mt1/2 b, m+1/2 _m+1/2 
x hi m+1/2 L ti "xx ) 


mia m+1/2 
bt"m,5 pa m+1/2 n gue » Sx ) 


m+1/2 


PT) A bet a 


tm+1/2 xxx 


La dimostrazione procede nello spirito delle prove dei Teoremi 3 e 4; il bel- 


lo è di controllare 


m+1/2 Mmt1/2 
Itm#1/20Amea/2 EI 


m+1/2 m+1/2 # Fade ) 


Ctraasali “ale 17 o 


US 


m+1/2 m+1/2 
Ct+1/2" dA] lp tie" 19 


HA 


b, m+1/2,2 10-b,2 
s 8* n+1/2! ex [+ Ct+1/2! A 


+ Ct hUA. 


peta/2 |? 10,2 
m+1/2!° o m+1/2 


+0 


Inoltre 


b, m+1/2 _m+t1/2 


t a (ny, » Sx 


+1/2 ) è maggiorato da 
m 


: db, m+1/2 44b 2 
$ t+1/2 h lesse IG A. 
Otteniamo 
2. |_mt1/2,2 d._mt1/2,2 
tr+:1/2%t1%m! + blî + Ctr+1/2" lei | 
2.,8 2b b+4 10-b 4 
s CA {h +tn+1/2! +tm+1/2" +tr+1/2! +(4t) } 
Con b = 10/3 


E 
i oa rita ab 
leggi le ol SCA CEIHTOg(zE )I(N°+(At) NT. 
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